Exact solution of double-5 function Bose gas through interacting anyon gas 
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Id Bose gas interacting through S, 8 and double-5 function potentials is shown to be equivalent 
to a S anyon gas allowing exact Bethe ansatz solution. In the noninteracting limit it describes an 
ideal gas with generalized exclusion statistics and solves some recent controversies. 
PACS numbers: 05.30.Jp, 03.70.+k 11.55.Ds, 71.10.Pm, 

The concept of particles with generalized exclusion statistics (GES) introduced by Haldane 0J has important 
consequences || in describing Id non Fermi- liquids ||, which in turn is believed to be related [Q to the edge excitations 
in fractional quantum Hall effect j5|. On the other hand, inspired by the success of the Chern-Simon theory, an attempt 
was made recently |6l to describe a Id ideal gas with GES in the framework of a gauge field model. However, in a 
subsequent paper [Q] the previous result was shown to be wrong and some other conclusion was offered. Our aim 
t-H ' here is to deal primarily with a Id Bose gas interacting through double-^ function potentials together with the well 
{-jTj. known 8 and derivative (5-function interactions. We show that this interacting model with several singular potentials 
is equivalent to a Id gas with GES (which we call anyon for brevity) interacting via 5-function potential only. This 
' 5-anyon gas is found to be exactly solvable by the coordinate Bethe ansatz (CBA) just like its bosonic counterpart, 
T^j- ' contradicting the common belief Q that the CBA is applicable only to models with symmetric or antisymmetric wave 
\ functions. Remarkably, at the limit of vanishing interaction the anyon gas becomes free and gauge equivalent to a 
related model proposed in [||. This shows that, though the explicit wave function and the A-body Hamiltonian of || 

j» are not exact, the conclusion it arrived at is basically correct. Therefore, while the error in the treatment of [|| was 

detected in [Q], the source of this error and the possible way to rectify it becomes evident from our result. 
We start with a Id Bose gas interacting through generalized 6 function potentials as 

T— I | AT 



H N = -^2dl k + 5{x k - xi){c + in(d Xk + d Xl )) ^ 5(^j - x k )S(xi - x k ) + j 2 ^ {S(x k ~ xi)) 2 . (1) 

k <k,l> <j,k,l> <k,l> 



^ 1 This model was briefly considered and readily discarded in Q as too difficult a problem to solve. Notice however 
Qh, that for 7 a = 0, a = 1,2, i.e. without the double-5 potentials, it has various exactly solvable limits. For example, 
for k = 0, c 7^ the model becomes the well known (5-Bosc gas [|l0|, while for n 0, c = it corresponds to Bose 
gas with S interaction p| . Both these cases are not only exactly solvable by CBA, but also represent quantum 
• i-H , integrable systems allowing i?-matrix solution. This can be proved through their connection with the quantum 
integrable nonlinear Schrodinger equation (NLSE) Jl^] and derivative NLSE [|l3]], respectively. Even the mixed case 
with k ^ 0, c ^ is solvable through CBA @Jlll|, though as a quantum model it does not allow a i?-matrix solution. 
Nevertheless for j a ^ 0, i.e. with the inclusion of highly singular double-(5 function interactions, the solvability of 
the model is completely lost and the application of the CBA becomes problematic due to the presence of tree-body 
interacting terms. We ask therefore, whether for some choice of the coupling constants j a other than zero, this 
difficulty could still be avoided and the solvability of the model be restored. We find the answer to be affirmative 
and in particular, for j a = k 2 the model becomes equivalent to a S function anyon gas, which appears to be exactly 
solvable similar to the well known bosonic case obtained at 7 a = k = 0. 

Instead of attacking model ([]]) directly, our strategy would be to transform it into some equivalent tractable 
problem. For this we notice that, parallel to the relation between the 5 Bose gas and the NLS model [Q, the 
generalized interacting bosonic model (|l]) can be considered to be the A-particle Hamiltonian of the nonlinear field 
model 
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H= I dx[: ((4V* + cp 2 + wsptyty- ~ = W(^p 2 ^)] , P = (2) 

involving bosonic operators ip^(y)] = S(x — y). In (^) we have chosen 71 = 72 = « 2 and introduced notation 

: : to indicate normal ordering (NO) in bosonic operators. Restricting now to the |TV > particle state and defining 
the TV-particle wave function as 

$(x n ,x i2 ,. . .,x iN ) =< 0\ip(x h )ip(x i3 ) • • ■ip(x iN )\N >, (3) 

we can generate all terms of (Q) starting from (|^). For example, the last term in (||): J dx^ p 2 ip) is equivalent 
to two normally ordered terms like J dx : (p 3 + p 2 J dy(S(x — y)) 2 ) : . The first one when acts on the state \N >, 
its three ip(x) operators in passing through the creation operators at points xj,Xk,xi in TV > would produce sum of 
terms with product of three 5 functions having arguments (x — Xj), (x — Xk) and (x — x{). On integration by x they 
would generate the double-i5 function potential 8{xj — Xf.)S(xi — Xk)- Note that this is a 3-body term and would not 
contribute in 2-body bosonic interactions. On the other hand, the second term acting on TV >, would give rise to the 
sum of terms like S(xk — xi)S(xk — Xk) ~ (S( x k — xi)) 2 - Similarly other terms with 6 and S function interactions are 
obtained in ([!]) from (0). 

Our next step is to define a gauge transformed operator 

4,(x) = e-' lK ^ 1{x ' mx ' )dx '^(x) (4) 

along with its conjugate tp^(x). We may check that the derivatives and products of the transformed operators are 
related to the old ones in the following way. 



i ; = : (^V) 2 : (5) 

: V>+V> x + iKiftp)(i) x - inpip) =: (^Va + Wp(tpH x - iptip)) : +k 2 (^ p 2 ip), (6) 

where : : stands for NO with respect to the transformed operator (|J), which not necessarily coincides with the 
bosonic NO as evident from (ph. Using these relations therefore one rewrites Hamiltonian (0) in the form 



H = dx 



(^ x + c(^W): (7) 



Note however that inspite of the same form as NLSE, (Q) is not same as the known model, since the fields involved 
are no longer bosonic operators but exhibit anyon like properties 

^{x 1 )^{x 2 ) = e iK <^- x ^^(x 2 )i'Hx 1 ), $( Xl )ft(x 2 ) = e- iK < x i-^ft(x2)i>{ Xl ) + S( Xl - x 2 ) (8) 

etc. where 

e(x — y) = ±1 for x > y, x < y and = for x = y, (9) 

which may be expressed also through the symmetrical unit-step function p5[ . This means that the bosonic commu- 
tation relation (CR) [ip(x), i/>'(y)] = S(x — y) remains valid at the coinciding points. These relations can be checked 
easily by using realization (|J) through bosonic fields. 

For finding TV-body Hamiltonian corresponding to (Q), we observe that operator ip(x) in passing through the string 

of anyonic creation operators in TV > would pick up first a phase e ' re 2_Ji<* e ^ x x ^ due to (@) and then leave a S(x — x k ) 
at Xk due to its standard CR at the coinciding points. The phase factor however is canceled subsequently when the 
associated i i/>' (x) also passes through the same creation operators and comes to the point Xk- This happens due to the 
opposite signs of the phases as seen from @. Therefore finally, similar to the bosonic model one obtains a 5 function 
interacting gas 

N 

H N = -J2d 2 Xk +c ]T S(x k - Xl ) (10) 

k <k,l> 



2 



However, in contrast to the standard case the wave function now exhibits a generalized symmetry 

$(n,...,a;i,...,a;j,...,afjv) =e l^*=*+i z^ =i+ i fc V$( Xl , . . . , Zj, . . . ,x i; . . . ,xjv) 



(11) 



dictated by the operator relations (j^) and the model is defined in an infinite-interval space R. It should be emphasized 
that due to the validity of bosonic CR at coinciding points in the anyonic relations (||), in the induced wave function 
(Ifi) the phase factor with e(xi — x^) vanishes at xi = x^ making it well defined. Note that the commutation relations 
(|) and the symmetry of the wave function ( |Ti"| ) for this Id model are remarkably consistent with the genuine anyonic 
behaviors in 2d |16|. Like the anyonic property, while ffl is q^ 1 — e* ree ( Xl ~ S2 )-symmetric, the wave function &(x\,X2) 
is g-symmetric, so that the 2 >-particle state remains invariant under permutation of coordinates. It is evident that 
the range of n is sufficient to be restricted in the interval of 2ir by choosing — 7r > k < ir and any nonzero value of 
k affects the symmetry under reflection of coordinates, close to the 2d anyon case. Another parallel of multi-anyon 
wave function Jn| can be observed in (p"T|), i.e. the phase factor appearing under exchange of two of its arguments 
also depends on the intermediate coordinates. 

We have converted thus the original eigenvalue problem related to (Q) to an equivalent one: Hn^(x%, . . . ,xn) = 
En<§>(xi, . . . ,xjy) for ( |io| ) acting on anyon-type wave function (|TT[). Since the Bethe ansatz solution is meant appar- 
ently for the (anti-)symmetric wave functions only ||, the present problem claims novelty. However, we find that the 
CBA is applicable here with equal success, if one modifies the Bethe ansatz for the wave function appropriately: 



$(xi, . . . ,x N ) = $a{xi, ■ ■ ■ ,x N )$ B (xi, . . .,x N ). 
Here $b is the symmetric function in the standard Bethe ansatz form p0| 



defined in the primary region R\ 

with e(xi — Xj) as defined in (^). 
boundary 



x\ < X2 < • ■ ■ < xn , while <&a is the additional anyonic part given as 

-i%A(x u ...,x N ) with 4 = ^^.^ 



i<j 



(12) 



(13) 



(14) 



Remarkably, the discontinuity in the derivative of the wave function ( p"2| ) at the 

(a,, - d ah )$\+ - (d xi - - c$| (15) 

with the notation \± =\ XI=X ± and |o = \ Xl =x k , determines the scattering amplitude in the same way as in case of 
(5-Bose gas. This becomes possible since using (fl2l ) , (|l4| ) along with (||) one obtains 





■ ,x N ) + 




• \ Xl=x + 


$(xi,.. 


.,X N ) |_ 




■■■' x n) \ Xl=x - 


Hx u . 


■ .,x N ) |o 


= e- i ^ s H B (x 1 ,... 


xn) \xi=x k 



(16) 

where S = A — e{xk — Xi). Contributions coming from the derivatives of other e(xi — Xj) factors (as S functions) in 
( |l5| ) clearly cancel each other, transforming it consequently to an equation only for the symmetric part of the wave 
function in the standard form 



(d Xl - d Xk )<f> B (x! 



,x N )\ + = c$ B (xi, ■ ■ ■ ,x N )\ , 



(17) 



but with modified coupling constant c = ic ^ s * . Note that the singularity of the original bosonic problem is reflected 
in the discontinuity of the anyonic wave function (|l2"|). Such discontinuity at the boundaries of different regions, 
though somewhat unusual in CBA, has been observed recently in other context fl7|| . In the present case, this also 
does not affect the physical picture, as evident from the reduction of the problem to ([17]) for continuous wave functions. 
Therefore following arguments of the bosonic model [ pl| , one can reduce ( p^ ) further to the region R\ involving only 
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adjacent k's, i.e. with x\ — > Xk+i and using the Bethe ansatz ( |l3| ) calculate the 2-particle scattering amplitude as 
gjflii+i _ ^' + !| = e~ j6l, +". Mark however, that in contrast to (To) the coupling constant in this anyonic case has 
been changed. At n — one recovers the bosonic case, while k — ±7r gives hard core repulsion c — ► oo, simulating 
fermionic model. In the present infinite-interval space the values of {k} have no restriction. However if we restrict to 
the interval < Xj < L, the boundary condition would be twisted like <&(x\+l, ■ ■ ■ , %n) = e lK ( N ~ 1 * > <&(xi, . . . , xn) 
and one would obtain the determining equations for k's as kj = — 4 X)fLi ®js + ^r n j + + 1 — j = 1, . . . , N 
with rij an integer. Here we have redefined = — Ke(j — s) = —0 S j to introduce an anyonic scattering phase, 
since such a particle in passing through others would pick up a phase e ±iK depending on its relative position. The 
energy eigenvalue of the system En = ^J> though has the same form as in bosonic model, acquires in effect 
different value due to changed coupling constant. We may mention here that the possibility of gauging away certain 
multi-species fermionic interaction by introducing a twisting in the boundary condition has been observed recently 
p8| . Similar spirit of the present result in a totally different context might therefore be an indication of a deeper 
generality. 

Now at c — > limit of ([!(]) we recover the results related to a GES ideal gas having properties like (||) and jll]). 
As we have shown, this free anyonic Hamiltonian would be equivalent to the ./V-bosonic model (0) at c = 0,7 a = n 2 . 
Establishing such a fact through direct gauge transformation of the wave function was the aim of Q. However this 
seems to be difficult to achieve, in which investigations were concentrated. We have avoided this difficulty by 
showing the equivalence through the gauge relation between the related field models (^) and (|^). Mark also that the 
field model (Eqn. (8)) of |j| differs from that (Eqn. (6)) of and ours by a crucial term. We note again that (@) 
even with c = involves 3-body terms, ignoring which naturally would make it not equivalent to the free anyonic 
Hamiltonian. In both Q and Jt| the authors worked with a 2-body bosonic Hamiltonian which could not give the 
right equivalence. 

The author likes to thank Partha Mitra for helpful discussions and Diptiman Sen for providing useful references. 
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